Spin Josephson effect with a single superconductor 
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A thin ferromagnetic layer on a bulk equal-spin-pairing triplet superconductor is shown to mediate 
a Josephson coupling between the spin f and 4- condensates of the superconductor. By deriving 
analytic expressions for the bound states at the triplet superconductor-ferromagnet interface, we 
show that this spin Josephson effect establishes an effective anisotropy axis in the ferromagnetic 
layer. The associated Josephson spin current is predicted to cause a measurable precession of the 
magnetization about the vector order parameter of the triplet superconductor. 
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Introduction. The complex relationship between super- 
conductivity and magnetism has motivated an enormous 
effort to understand the properties of heterostructure in- 
terfaces between ferromagnets (FMs) and spin-singlet su- 
perconductors (SSCs)ji A remarkable feature of such de- 
vices is the existence of proximity-induced spin-triplet su- 
perconducting correlations due to the exchange splitting 
in the FM, which are responsible for the anomalous dy- 
namics of the barrier magnetization in an SSC-FM-SSC 
Josephson junction^ Due to the intimate connection be- 
tween ferromagnctism and triplet superconductivity, it 
is natural to consider what results if the SSC were re- 
placed by a triplet superconductor (TSC). This question 
is of fundamental interest, as the intrinsic spin struc- 
ture of the Cooper pairs in a TSC allows us to anticipate 
an unconventional and unique interplay with magnetism, 
which may be unambiguous signatures of the triplet pair- 
ing state in proposed TSCs such as LiFeAs£ For exam- 
ple, bulk spin supcrcurrents are known to be possible in 
TSCs^r— and some proposals for their realization require 
FM elements. The response of the FM component of the 
device to the spin supercurrent, however, has yet to be 
investigated. The recent fabrication of superconducting 
thin films of the suspected TSC Sr2Ru04 is an important 
step towards the creation of TSC-FM heterostructures^ 
and so a deeper understanding of the physics of TSC-FM 
interfaces is timely. 

In this paper we show that a thin FM layer on a bulk 
cqual-spin-pairing TSC produces a spin Josephson ef- 
fect by coupling the spin f and J, Cooper pair conden- 
sates. The physical mechanism is the spin-dependent 
phase shift acquired by a Cooper pair undergoing spin- 
flip reflection at the FM interface, which acts analogously 
to the phase difference in a Josephson junction. Making 
only the assumption of spatially constant order parame- 
ters, we solve the Bogoliubov-de Gennes (BdG) equations 
for the bound states at the TSC-FM interface. Using this 
we calculate the free energy of the interface, revealing 
that the spin Josephson effect creates an effective hard 
or soft axis within the FM layer, depending upon the 
orbital structure of the TSC gap. Finally, we obtain a 
general expression for the spin current, thereby showing 
that it exerts a measurable torque on the FM moment. 
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FIG. 1. (Color online) Schematic diagram of the device stud- 
ied here. The three different choices for the TSC orbital are 
shown: for the p z and p y cases, the white and black lobes 
indicate opposite signs; the arrow in the p z + ip y shows the 
direction of increasing phase. 



We propose this effect as a test for triplet pairing. 

Model system. We consider a FM layer of width L on 
a bulk TSC, separated by an atomically thin insulating 
layer, see Fig.Q] The BdG equation for the quasiparticle 
states with energy E is 
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where the caret indicates a 2 x 2 matrix in spin space. 
The wavefunction ^(r) is only non-zero for z < L. The 
non-interacting Hamiltonian Hq(t) is 
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Here \xb is the Bohr magneton, and we make the simpli- 
fying but nonessential assumptions that the radius of the 
Fermi surface &f and the effective mass m in the normal 
state of the TSC and the FM are the same. The insulat- 
ing layer is modeled as a 5-function of strength U. The 
last term in Eq. ([2]) is the energy due to the exchange 
field H cx in the FM. For an incompletely polarized FM, 
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the magnetization m = |m|[cos(o;)e a ; + sin (a) cos(?7)ej, + 
sin(a) sin(?7)e 2 ] is related to H cx by Luttinger's theo- 
rem. We assume a two-dimensional system, so that the 
majority-spin (parallel to H ex , s = +) and minority- 
spin (antiparallel to H ox , s = — ) Fermi surfaces have 
radius kp , 



y(T+ s\)kp, where A = /is|H 0X |//i < 1. 
This gives |m| = y? B m\H. 0X \/ i:h 2 . For a half-metallic FM 
there is a single Fermi surface of radius V2kp, and the 
exchange splitting is fixed by the details of the system. 

The gap matrix is A(r) ~ 0(— z)i[a ■ d]a y where 
d = Ad is the vector order parameter, assumed con- 
stant throughout the TSC. A is an operator, which for 
Cooper pairs in a relative p-wave orbital state has the 
real space form A = — iA(T)n • V/fc^. The gap mag- 
nitude A(T) is assumed to have weak-coupling temper- 
ature dependence. The unit vector n defines the orbital 
state: n = e 2 for p 2 -wave; n = e y for p^-wave; and 
n = e 2 + ie y for (p z + ip y )-w&ve. The last choice is of 
greatest relevance to Sr2Ru04 and LiFeAs^ while the 
others have been proposed for (TMTSF) 2 X (X = PF 6 , 
C104)J£ In the following it is convenient to express the 
gap in terms of the Fourier transform of A, which is writ- 
ten Ak = A(T)n ■ k/kp. We fix d = e x which defines a 
TSC where the Cooper pairs have z-component of spin 
S z = ±h but the condensed part of the system is unpo- 
larized. Below we show that only the angle between d 
and m is relevant for the spin Josephson effect, and so 
other orientations of d do not result in new physics. 

Bound states. We seek solutions of Eq. (Q]) for states 
bound to the FM layer. The wavefunction of such a state 
has the general form , 5(kj|;r) = ^E , TSc(k|| ; r)0(— z) + 
5 , FM(k|| ; r)6(z) and satisfies lim^-oo \I/(k||;r) = and 
^(kn ; r)| 2= L = 0. The momentum component parallel 
to the interface, k||, is a good quantum number due to 
translational invariance. Solving the Andreev equations 
in the TSC ) 11 ' 12 we make the ansatz 

*TSc(k||;r) = [ai,a*a(ki;r) + a 2 ,a* CT (k 2 ;r)] , 

<T=t4 

(3) 

where the spinors are given by 

* t (k;r) = (l, 0, 7 (k), 0) T e ikr e^ , (4) 

* 4 .(k;r)=(0 > 1, 0, - 7 (k)) T e^'e"** , (5) 

with 7 (k) = -[g + tsgn (fc 2 )V|A k | 2 -£ 2 ]/A k and K k = 
{rn/h 2 \k z \)\J\A\ s \ 2 — E 2 . The wavevectors appearing 
in Eq. (j3|) are defined by ki = (k||,fc 2 ), k 2 = (k||, — k z ). 
Note that |A kl | = |Ak 2 | = |A k|| | for the orbital symme- 
tries considered here. 

Depending upon the value of k|| , we have either prop- 
agating or evanescent solutions in the FM layer. In the 
case when there are propagating solutions in both spin 
channels we have 

*FM(k||;r) = J2 {^sin(fc e , s [L-z])e ik n- r $e, s 

s=± 

+ b htsS in(k h<s [L-z])e ik r^ htS } (6) 



where the electron and hole spinors are defined by 

$e, S = (w B , x s , 0, 0) T , $ h , s = (0, 0, w* s , x s f , (7) 



with w s = s(cosa — isinacosrj)/^/l — s sin a sin rj, x s = 
i/l — ssinasin?7, and the wavevector fc e (fc),s for electrons 
(holes) is fc e(h)iS = [fc|(l + sX) - |k|| | 2 + (-)2m£/fi 2 ] 1 / 2 . 
If the radicand is negative, only evanescent solutions are 
possible; in this case we replace k e (h),s - > * K e(/i),s where 
K e(h),s is the inverse decay length. 

The coefficients in Eq.s © and © are chosen so that 
at the TSC-FM interface the wavefunction is continuous 
^(kn ; r)| 2=0 - = ^(ky ; r)| 2=0 + , and its derivative obeys 
9^(k||;r)| z=0+ - 0,¥(k,|;r)|, =o - = 2Z¥(k,|;r)|, =0+ 
where Z = rnU/h 2 . The values of E for which the 
determinant of the resulting system of equations van- 
ishes define the bound-state energies. Explicit expres- 
sions for the bound-state energies can be found when 
the ^-dependence of the wavevectors is neglected, i.e. 
k e ,s ~ fcft,s ~ k s . This approximation is valid for a thin 
FM layer'such that (k e , 8 - k h . s )L w 2EL/hv F . s < 
For a weakly to moderately polarized FM layer we have 
Fermi velocities Vp t + ~ vp t - ~ 10 ms~ , and so for 
E < max{|Ak|} ~ O.lmcV (T c ~ IK) we require thin 
layers less than about 100 unit cells thick. In this limit 
we obtain the non-degenerate bound states: 
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Here we have 
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where A = A/^k 2 F - |k||| 2 (A = k s ,Z). 

The bound-state energies Eq. (|5J) are a central result of 
our paper. They originate due to multiple Andreev reflec- 
tions within the thin FM layer, which Josephson-couplc 
the S z = ±h condensates in the TSC. The same physical 
mechanism is responsible for the formation of Andreev 
bound states (ABSs) at the tunneling barrier in a Joseph- 
son junction ! 12 ' 14 ' 15 Remarkably, the states Eq. (JSJ arc 
identical to the (spin degenerate) ABSs in a short Joseph- 
son junction of transparency between p 2 -wave TSCs 
with phase difference Acj) = 2a [Eq. (|8a|) ] . between p y - 
wave TSCs with Ac/) = n + 2a [Eq. ([8b| ]. and between a 



(p z + ip y )-wave and (p z — ip y )-~wave TSC with Acf> = 2a 
[Eq. ([5c])]. Since the form of the bound-state energies is 
fixed by the bulk pairing symmetry, our results should 
be robust to a self-consistent calculation of the gapJ^ 

The spin Josephson coupling can also be understood 
at a more fundamental level: in the BdG Hamiltonian 
the S z — ±h condensates in the TSC are independent of 
one another. Tunneling of a Cooper pair between the two 
condensates is made possible by the FM layer, where the 
coupling to the FM moment allows an incident Cooper 
pair with spin ah to be reflected with spin —ah. As a 
result of this process, the Cooper pair acquires a phase 
shift A0kn + 7T — 2a a. The last terms are due to the 
spin flip itself, and are primarily responsible for driving 
the spin current. A#kii = ar g{Ak 2 } — arg{Ak x } is the 
phase shift due to the orbital structure of the TSC: for 
p z -wave orbitals we have A^ky = tt for all k|| ; in the 
Py-w&ve case we have A(9k|| = 0; and the superposition 
of these two orbitals in the (p z + ip y )-w&ve TSC gives 

A#k|| = t — 2arccos(yl — |kjj \ 2 /k 2 F ). A6*k y accounts for 
the 7T phase difference between the bound states in the 
p z and p y cases, and the apparent sign reversal of the p y 
component in the p z + ip y bound states. 

We note that unlike the standard Josephson effect, 
where the phase difference between the two superconduc- 
tors drives the supercurrent, the spin Josephson effect is 
due to phase shifts picked up during the tunneling pro- 
cess itself. The phase difference between the spin-up and 
spin-down condensates is fixed by the orientation of d, 
which is unaffected by the FM layer. 

Free energy. The assumption of spatially constant or- 
der parameters allows us to write the free energy due to 
the spin Josephson coupling in the TSC-FM device in 
terms of the bound states^ 



F = -\k B T ^ ^ In [2 cosh(/3£„, k|| /2)] + F . 



(10) 



n— ± kn 



where Fq is independent of a and includes the interaction 
energy in both the TSC and the FM, as well as the con- 
tribution from continuum states. We plot the free energy 
difference AF(a) = F(a) — F(a = 0) per interface-unit- 
cell area a x a y as a function of a in Fig. Ufa). As can 
be seen, F takes a minimum as a function of the angle 
a: regarding d as a fixed property of the bulk TSC, the 
spin Josephson coupling therefore establishes a preferred 
orientation for the magnetization of the FM layer. There 
is a direct analogy to a short Josephson junction, where 
the free energy due to the Josephson coupling takes a 
minimum as a function of Arj)^ 

In the p z -w&ve (p y -w&ve) case, the formally equivalent 
Josephson junction with A(f> = 2a (A</> = n — 2a) has 
time-reversal symmetry and so the Josephson free en- 
ergy is minimized at A<j) ~ 0. This implies that in the 
TSC-FM device the free energy always has a minimum 
at a = (a = %/ '2). The broken time-reversal symme- 
try in the {p z + ip y )-wave case, however, means that the 
stable value of a is determined by the details of . 




FIG. 2. (Color online) (a) The free energy difference AF(a) 
and (b) the z-component of the spin current js,z per interface 
unit cell area as a function of a for the three choices of orbital 
wavefunction. We take A = 0.05, T = 0.4T C , n = 0, L = 10a z 
and Z = 1. Ao is the T — gap magnitude. 



Specifically, the p z and p y components of the gap favor 
minima at different values of a: if is peaked near 
|k|j| = kp (0), the p y -componcnt (p z -component) domi- 
nates and the configuration with a = 7r/2 (0) is stable; 
for more complicated , the competition between the 
gap components may stabilize the system at a ^ 0, 7r/2. 
For weak magnetization strengths, the free energy mini- 
mum is usually located at a = 0. 

From Fig. [2ja) we see that to excellent approximation 
F oc cos 2a. In writing an effective free energy for the FM 
layer, we can therefore account for the spin Josephson 
effect by including a term Fj = f s (d ■ m) 2 , i.e. for f s < 
(f s > 0), d defines an effective easy (hard) axis in the FM 
layer. Since the sign of f s is determined by the orbital 
state of the TSC, this reveals a novel type of spin-orbit 
coupling between the TSC and the FM. 

Spin current and magnetization dynamics. The zero- 
bias charge current Ij in a Josephson junction is given by 
Ij = (2e/h)dF/dAcf). We now show that in our device 
there is a spontaneous spin current which can be similarly 
expressed derivative of F with respect to a. 

Our starting point is the continuity equation for the 
spin 



h d _ K OF 
-M = M x — — 

qHB dt dM 



(11) 



where g is the gyromagnetic ratio and M is the total 
moment of the FM layer. The vector notation for the 
spin current J s refers only to the polarization; the direc- 
tion of the spin current is normal to the interface. We 
omit gradient terms in Eq. (|lip as the FM layer is con- 
sidered to be thin compared to the coherence length of 
the Cooper pairs. The free energy in our problem has 
the form F = F(|M|, d • M), allowing us to write 



= M x 

= M x 

_ .dF 
^ da 



d\M\ dF ^ dd ■ M OF 
dM d\M\ ^ 

dF 



dM dd-M 



1 . ... dF 

M m + w\ [d ^ d - M)M] d^Ma) 



(12) 
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where p = d x M/|M| sin a is a unit vector which points 
in the same direction for all M lying in a fixed plane 
containing d. Inserting Eq. (fit)]) into Eq. (|12p we obtain 

1 dE 

J.s = -P i E E ^M/^n.k,, 12) . (13) 

n=± ky 

This closely resembles the Beenakker-van Houten for- 
mula for the charge current in a short Josephson junc- 
tion.— It reveals that the spin current in our device is 
due entirely to resonant tunneling between the two spin 
condensates through the bound states Eq. ©. Equa- 
tion ([13]) gives identical results to the Furusaki-Tsukada 
technique,— which expresses the spin current in terms 
of the Andreev reflection coefficients. We show the spin 
current as a function of a in Fig. [2][b). 

If the magnetization is prepared with < a < 
7r/2, Eq. (|lip predicts that the spin current will exert 
a torque on M, causing it to precess about d. Writ- 
ing M = AL/iBpnm and J s = Aj Si where A is the 
area of the TSC-FM interface, and p = (n+ — n_)/n 
is the polarization of the FM, we find the precession fre- 
quency to be il j = 2g cos(a)max{|j s |}/7mpL. To esti- 
mate f2j, we assume a weakly polarized FM, A = 0.05, 
with n = 1 electron per unit volume v = a x a y a z , 
max{|j s |} = 0.025A /(a x a y ) [see Fig.^b)], L = 10a z , 
and T c = IK. We hence find ftj = 15cos(a) GHz, which 



is measurable by ferromagnetic resonance (FMR) exper- 
iments. As the spin Josephson effect cannot occur for 
a SSC, the observation of this precession would be very 
strong evidence of a triplet pairing state, although the 
precession effects due to multiple FM domains would 
have to be ruled outJi Similarly, it is also necessary 
to examine the effect of chiral domains of the gap in 
the candidate material Sr2Ru04>i£ Gilbert damping and 
anisotropy effects in the FM layer must also be included 
in a complete description of the magnetization dynamics, 
but do not change the derivation of Eq. (fT3"f . 

Conclusions. In this paper we have demonstrated that 
the spin structure of the Cooper pairs in a TSC permits 
the occurrence of a spin Josephson effect without the 
need for a second superconductor. We have proposed 
that a thin FM layer on a bulk TSC can realize this ef- 
fect. In turn, the spin Josephson coupling establishes 
an effective easy or hard axis in the FM layer, depending 
upon the orbital symmetry of the TSC gap. Furthermore, 
the Josephson spin current causes the magnetization to 
precess about the d vector with a frequency that is acces- 
sible to FMR, realizing a possible experimental signature 
of the triplet state. 
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